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We evaluated rates of transversing muons, muon-induced fast neutrons, and production of 40 C1 
and other cosmogenically produced nuclei that pose as potential sources of background to the physics 
program proposed for an argon-based long baseline neutrino experiment at the Sanford Underground 
Research Facility (SURF). The Geant4 simulations were carried out with muons and muon-induced 
neutrons for both 800 ft (0.712 km.w.e.) and 4850 ft levels (4.3 km.w.e.). We developed analytic 
models to independently calculate the 40 C1 production using the measured muon fluxes at different 
levels of the Homestake mine. The muon induced 40 C1 production rates through stopped muon 
capture and the muon-induced neutrons and protons via (n,p) and (p,n) reactions were evaluated. 
We find that the Monte Carlo simulated production rates of 40 CI agree well with the predictions from 
analytic models. A depth-dependent parametrization was developed and benchmarked to the direct 
analytic models. We conclude that the muon-induced processes will result in large backgrounds to 
the physics proposed for an argon-based long baseline neutrino experiment at a depth of less than 
4.0 km.w.e. 

PACS numbers: 13.85.Hd, 23.40.-s, 25.40.Fq 



I. INTRODUCTION 

Experiments for the past several decades require mod- 
ification of the Standard Model to incorporate the unex- 
pected neutrino properties and fundamental characteris- 
tics [H-Q • For instance, neutrino flavor mixing was found 
to be responsible for the phenomenon of neutrino oscilla- 
tion lS|] occurring between three generations in which 
a complex phase (Sep) signifies the amount of violation 
of the charge-parity (CP) symmetry, which is unknown. 
The value of the mixing angle, #13, between the first gen- 
eration and the third generation, is also unknown, but an 
upper limit of <10° is given by CHOOZ experiment j| 
at the 90% C.L. The sign of the mass difference, Ato 2 3 , 
which represents the ordering of the mass eigenstates, 
remains unknown as well. These three unknown param- 
eters are intended to be addressed with the next genera- 
tion of neutrino oscillation experiments. 

Recent studies of neutrino properties using neutrino 
beams have demonstrated that good sensitivity to CP 
violation and the mass hierarchy can be achieved by mea- 
suring v e appearance using a very long baseline beam 
with massive detectors [l0| assuming the value of 6*13 > 
1°. The recent values of 6*13 measured by MINOS and 
T2K shows a non-zero value [ll|> [l2j that is particularly 
interesting to the measurements of CP violation and mass 
hierarchy. The proof of CP violation in the lepton sector 
and the knowledge of the value of Sep are crucial to un- 
derstanding the origin of the baryon asymmetry in the 
universe, providing a strong implication of leptogenesis 
that is responsible for the observed baryon asymmetry of 
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the universe [13| . Simultaneously, the neutrino mass hier- 
archy is of great im por tance for neutrinoless double-beta 
decay experiments [14| and could shed light on possible 
flavor symmetries. 

When measuring the value of #13, CP phase, Sep, and 
the neutrino mass hierarchy with conventional neutrino 
beams, a key process of new discovery in neutrino os- 
cillation is v e appearance. With appropriate de- 
tector design and adequate control of environmental fac- 
tors, specfically proper shielding, the long baseline neu- 
trino experiment (LBNE) [IH, [l6| is capable of support- 
ing an extremely rich program of high energy physics 
and particle-astrophysics including proton decay, astro- 
nomical neutrinos, and tests of fundamental physics and 
the Standard Model. The two detector technologies be- 
ing considered are: 1) an active finely grained liquid ar- 
gon time-projection-chamber (LAr-TPC) and 2) a water 
Cerenkov detector. Both detector technologies can sup- 
port this wide range of physics goals probing the Stan- 
dard Model and searching for physics beyond the cur- 
rent models. However, muon-induced background may 
constrain the sensitivity of the proposed experiment, in 
particular, the proposed galactic and relic supernova neu- 
trino detection with liquid argon detector. Though the 
background can be measured for beam neutrino physics 
with the beam-on or -off, the fluctuation of background 
events can be problematic if these events occur on the 
same order of magnitude as anticipated physics signal. 
In addition to the beam contamination backgrounds, the 
main sources are the muon-induced processes. In this 
paper, we report the study of the muon-induced back- 
grounds for an argon-based detector. 

Muons and muon-induced fast neutrons entering the 
detector from the surrounding rock can result in trouble- 
some experimental backgrounds. While through- going 
muon events in the detector have energy depositions and 
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can be easily identified, muons that traverse the rock 
near the detector, or traverse very small distances in the 
detector ("corner-clippers"), can result in limiting back- 
grounds. Moreover, background events can be produced 
by the muon-induced fast neutrons entering the detec- 
tor while the parent muons pass completely through the 
surrounding in-active materials. The muon-induced neu- 
trons can undergo (n,p) reaction with Ar, and negative 
muons capture on 40 Ar. Both reactions create unwanted 
40 CI, which can be a background for the following re- 
actions: z7 e + 40 Ar ->• e-+ 40 Cl, zy e + 40 Ar^ e -+ 40 K, v x + 
40 Ar — > v x + 40 Ar*, and v Xt x + e - — » v x + e~, where x = 
e, jit, t. This is because the decay Q value of 40 CI is 7.48 
MeV which is above the proposed detection threshold of 
5 MeV or 6 MeV, and the half-life of 40 CI is 1.35 min- 
utes making it difficult to correlate with muons. There 
can also be additional radioactive isotopes produced by 
muon-induced processes in the argon target. Those ra- 
dioactive isotopes can be part of the background for the 
proposed physics channels. Therefore, it is necessary to 
evaluate the muon-induced processes and the production 
rate of the cosmogenics in the detector volume. 

Due to the lack of direct detection measurements of the 
muons and muon-induced products and parameters for a 
given depth with a given detector, such an evaluation has 
not been modelled in full. In this paper, we present sev- 
eral parametrization functions that estimate the muon- 
induced fast neutron energy spectrum and the stopping 
muons as a function of depth. Using these parametriza- 
tion functions, we simulated cosmogenic production rates 
in the proposed LBNE LAr detector for two depths with a 
well known Geant4 package [HI, Geant4.9.5 with shield- 
ing physics list. The simulated results were compared to 
the predictions by the developed analytic models. Since 
the muon-induced processes are strongly depth depen- 
dent [ItJ , we establish a depth-sensitivity relation for an 
argon-based detector by calculating the production rates 
of muon-induced processes as a function of depth. 



II. EVALUATION OF MUON-INDUCED 
BACKGROUND 

At sufficiently high energies, radiative processes be- 
come predominant in energy loss for muons. Many sub- 
sequential particles can be produced by muon-induced 
radiative processes. Cosmogenic radioactive isotopes can 
be produced by muons and muon-induced sequential par- 
ticles including neutrons, protons, pions, gamma rays, 
etc. Two cases are considered in the following evaluation: 
1) muons transversing the detector and 2) muon-induced 
neutrons entering the detector from the experimental 
hall. The former process creates 40 C1 and other radioac- 
tive isotopes through negative muon capture, (n,p) reac- 
tion, and (p,n) reaction, etc. The latter produces back- 
grounds mainly through high energy events, ttq and 40 C1 
productions. In order to understand the production rates 
and their corresponding mechanisms, we performed a full 



Geant4 Monte Carlo simulation and developed analytic 
models. We elaborate on the evaluation processes in the 
following subsections. 



A. Muon-induced Background from the Geant4 
Simulation 

High energy cosmic-ray muons can penetrate rock over- 
burden to reach an underground laboratory with the 
surviving muons generating neutrons in the surrounding 
rock. Those neutrons are unwanted particles that could 
produce background events for low-background experi- 
ments searching for rare event physics. The intensities 
of the residual muons and the muon-induced neutrons 
depend strongly on the depth of the underground detec- 
tor. This is particularly important for detecting super- 
nova neutrinos that contribute to a signal range of a few 
MeV to a few tens of MeV, which can often be dom- 
inated by the muon-induced backgrounds depending on 
the depth. Therefore, the depth-sensitivity relation needs 
to be understood in order to choose an appropriate depth 
at the Sanford Underground Research Facility (SURF) to 
estimate the experimental limitations imposed by these 
backgrounds and to optimize the detector design for the 
full range of possible physics programs. Because both the 
800-ft level and the 4850-ft level are being considered for 
a far detector, we have conducted Geant4 Monte Carlo 
simulations to understand the cosmogenic production in 
the detector at both levels. 

In the Geant4 simulation, a simple geometry with a di- 
mension of liquid argon, 45.6 m (width) x 22.4 m (length) 
x 14.0 m (height) [lfl, was placed in a stainless steel con- 
tainer with a thickness of 1 cm. Both muons and neutrons 
were generated on a very thin sheet (100 m (width) x 50 
m (length)) right above the upper surface of the stain- 
less steel tank. Since the thin sheet is two times larger 
than the size of the detector in length and width, both 
muons and neutrons generated in the sheet can enter the 
detector with angular distribution of sec(9). There is a 
wide range of thetas that are enabled by taking a large 
production sheet and better represent the range of an- 
ticipated underground muons. Note that we assume the 
parent muons and daughter neutrons have the same an- 
gular distribution. The residual muon energy spectrum 
was obtained from equation (8) of Ref. [17[ . The normal- 
ization constant was determined using equation (4) in 
Ref. [l7| . The muon-induced neutron energy spectrum 
was given in equation (14) in Ref. [ItJ ■ However, this 
equation cannot be directly used to generate neutrons 
for a given depth without knowing the associated param- 
eters. We adopt the following procedures to obtain the 
neutron energy spectrum as a function of depth. First, 
the average neutron energy as a function of depth was 
studied using the simulated results for various depths in 
Tables I and VII from Ref. [l?} and a measured surface 
data point from Ref. [lj| . Fig. Q] shows the fitted curve. 
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FIG. 1: Average neutron energy as a function of depth. 

The fitted function can be expressed as 

<E n >= 113 ■ [1 - exp(-0.545 ■ h)], (1) 

where h is the depth in km.w.e. At 800 ft, < E n > 
is 36.3 MeV, and at 4850 ft, < E n > is 102.2 MeV. It 
is noticed that the average neutron energy at the depth 
of Gran Sasso is about 91 MeV [l7| ■ This is not much 
different from the average neutron energy predicted for 
the 4850-ft level at Homestake. This is because the pro- 
duction mechanisms are primarily sensitive to the over- 
burden of rock and do not depend much on site specific 
details. Therefore, we then used the parameters in Table 
VII of Ref . [13| to generate neutrons for the depth of Gran 
Sasso. The normalization is done with equation (13) in 
Ref. [17| that predicts the muon-induced neutron flux as 
a function of depth for depths greater than 1.6 km.w.e. 
One caveat in this method is the ignorance of rock com- 
position, which could make a difference of up to about 
35% [17] . The muon-induced neutron energy spectrum 
at the 800-ft level was obtained using a scaling method. 
We scaled the neutron energy spectrum from the depth 
of 4850 ft to the depth of the 800 ft using a scaling fac- 

tOT ' (!J^) Q fe^> where < ^,800 >, $,,800, and 
< -2^,4850 >, $,,4850: are the average muon energies and 
the total muon fluxes for the 800-ft level and the 4850-ft 
level at Homestake, respectively, a — 0.73 is a constant. 
Note that the primary concerns of the cosmogenic pro- 
duction are 40 CI and 40 K through negative muon capture, 
(n,p), and (p,n) reactions. Since the reaction threshold 
of 40 Ar(n,p) 40 Cl requires neutrons with kinetic energy 
greater than 6.87 MeV, this cuts the majority of neu- 
trons induced by natural radioactivity in rock through 
(a,n) reactions. Therefore, we neglected the calculation 
of 40 C1 production by (a,n) neutrons. 

1. Muon-induced Background at the 800-ft Level 

The muons that survived the 800-ft rock at Homes- 
take have an average energy of ^97 GeV and a flux of 



6.3xl0 -6 cm -2 s~ 1 from equation (9) and equation (4) in 
Ref. [13 • These high energy muons passing through the 
surrounding rock of a laboratory will generate fast neu- 
trons. The emerging neutrons in an experimental hall 
have an average energy of ~36 MeV with a total flux 
of 3.2xl0 _ cm, s . The neutron energy spectrum ob- 
tained from the methods described above is shown in 
Fig. [5] We simulated the muons passing through the de- 
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FIG. 2: Local fast neutron energy spectrum at the 800-ft 
level used in the Monte Carlo simulation. 

tector and the muon-induced neutrons entering the de- 
tector from rock. Fig. [3] shows a visible energy spectrum 
induced by the muon-induced neutrons entering the de- 
tector from the surrounding rock. The event rate above 5 
MeV in a 20 kton detector is ~0.28 Hz. The direct muon 
rate is ^88 Hz. These rates are high enough to poten- 
tially swamp any signals from a galactic supernova (~44 
Hz estimated using Ref. [2(|) or even from a neutrino 
beam (^75 events per year for u e appearance estimated 
with Ref. (U). 




Energy Deposition (MeV) 

FIG. 3: Fast neutrons induced visible energy spectrum in 
the detector at the 800-ft level. A fiat ionization efficiency of 
25% was applied to the nuclear recoil events in this plot. 

We summarize the 40 CI production rates from the 
above two sources in Table Q] It is also interesting to 
show the overall cosmogenic production in the detector. 
As shown in Table [TTI and Fig.[U the cosmogenic isotopes 
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TABLE I: 40 C1 production rates in the detector (20 kton) at 
the 800-ft level. 



From n simulation 


From n simulation 


Produced by 


Rate per day 


Produced by 


Rate per day 


Muon Capture 


27344 


Secondary fi 


45 


Secondary n 


40587 


Neutrons 


3667 


Pions 


249 


Pions 


1.4 


Others 


83 


Others 


< 1 


Total 


68163 


Total 


3714 



induced by muons range from P to Ca, a total rate of 19 
Hz. Some of them can be background to the proposed 
physics channels. 



TABLE II: Additional significant 
rates in the detector (20 kton) at the 



cosmogenic production 
800-ft level. 



Isotope 


Produced by 


Rate per day 


Q (MeV) 


tl/2 


30 p 


Spallation 


9020 


4.23 


2.5 m 


32p 


Spallation 


20900 


1.71 


14. 3 d 


33p 


Spallation 


30100 


0.25 


25.3 d 


34p 


Spallation 


12090 


5.4 


12.4 s 


35p 


Spallation 


7500 


4.0 


47. 2 s 


36p 


Spallation 


1190 


10.4 


5.6 s 


37p 


Spallation 


550 


7.9 


2.3 s 


31g 


Spallation 


5500 


5.4 


2.6 s 


35 g 


Spallation 


215500 


0.17 


87.5s 


37 g 


(n,a) 


31500 


4.9 


5.1 m 


38 g 


Spallation 


11500 


2.9 


170 m 


39 g 


Spallation 


850 


6.6 


11.5 s 


33 CI 


Spallation 


670 


5.6 


2.5 s 


34 C1 


Spallation 


8700 


5.6 


32 m 


36 CI 


Spallation 


1005000 


0.7 


3.1xl0 5 y 


38 C1 


Spallation 


110000 


4.9 


37.24 m 


35 Ar 


(n,6n) 


7100 


6.0 


1.8 s 


37 Ar 


(n,4n') 
(n,2n') 


21000 


0.8 


35 d 


39 Ar 


91000 


0.57 


269 y 


41 Ar 


capture 


45100 


2.5 


109 m 


38 K 


Spallation 


650 


5.9 


7.6 m 


40 R 


(P> n ) 


6500 


1.3 


1.28 x10 s y 


Total 




1641920 
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FIG. 4: Muon-induced cosmogenic production at the 800-ft 
level. 



than 6.87 MeV is about 1.6x 10 _10 cm _2 s _1 




2. Cosmogenic Production at the 4850-ft Level 
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FIG. 5: Local fast neutron energy spectrum used in the 
Monte Carlo simulation. 

Similar to the simulation for the 800-ft level, we simu- 
lated both the residual muons crossing the detector and 
the muon-induced neutrons entering the detector from 
rock. The visible energy spectrum induced by the muon- 
induced neutrons entering the detector from the sur- 
rounding rock is shown in Fig. [6j The event rate above 5 
MeV is ~0.001 Hz. The direct muon rate is -0.05 Hz. 

Table Mil shows the simulated 40 CI production rates 
from two sources. 



The residual muons at the 4850-ft level have an av- 
erage of —321 GeV using the equation (9) in Ref. fl7j . 
The total flux is predicted to be 4.4x 10 _9 cm _2 s _1 in 
Ref. |17| . The neutrons that are produced in the rock by 
these residual high energy muons entering the experimen- 
tal hall were simulated by Mei & Hime in great detail [llj ■ 
We obtained a muon-induced neutron energy spectrum 
using the method described in Section II A. Fig. [5] shows 
the neutron energy spectrum. The total flux at the 4850- 
ft level is 5.4xl0~ 10 cm~ 2 s _1 . The production of 40 C1 is 
via (n,p) reaction on 40 Ar requiring a threshold of 6.87 
MeV. The total neutron flux with neutron energy greater 



TABLE III: 40 CI production rates in the detector (20 kton) 
at the 4850-ft level. 



From simulation 


From n simulation 


Produced by 


Rate per day 


Produced by 


Rate per day 


Muon Capture 


17.5 


Secondary fi 


0.43 


Secondary n 


54.4 


Neutrons 


9.3 


Pions 


0.33 


Pions 


0.016 


Others 


0.04 


Others 


0.002 


Total 


72.3 


Total 


8.41 



Similar to the above discussion of the cosmogenic pro- 
duction for the 800-ft level, we show overall cosmogenic 
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FIG. 6: Fast neutrons induced visible energy spectrum in 
the detector at the 4850-ft level. A flat ionization efficiency 
of 25% was applied to the nuclear recoil events in this plot. 
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FIG. 7: Muon-induced cosmogenic production at the 4850-ft 
level. 



isotopes induced by muons for the 4850-ft level in Ta- 
ble |IV] and Fig. [7] Note that the production rate is rela- 
tively small at this depth. 



TABLE IV: Additional significant 
rates in the detector (20 kton) at the 



cosmogenic production 
4850-ft level. 



Isotope 


Produced by 


Rate per day 


Q (MeV) 


tl/2 


30 p 


Spallation 


9.6 


4.23 


2.5 m 


32p 


Spallation 


22.2 


1.71 


14. 3 d 


33p 


Spallation 


31.9 


0.25 


25.3 d 


34p 


Spallation 


12.8 


5.4 


12.4 s 


35 p 


Spallation 


8.0 


4.0 


47. 2 s 


36 p 


Spallation 


1.3 


10.4 


5.6 s 


37p 


Spallation 


0.6 


7.9 


2.3 s 


31g 


Spallation 


5.8 


5.4 


2.6 s 


35 g 


Spallation 


228.5 


0.17 


87.5s 


37g 


(n,a) 


33.4 


4.9 


5.1 m 


38 S 


Spallation 


12.2 


2.9 


170 m 


39g 


Spallation 


0.9 


6.6 


11.5 s 


33 CI 


Spallation 


0.7 


5.6 


2.5 s 


34 C1 


Spallation 


9.2 


5.6 


32 m 


36 CI 


Spallation 


1065.7 


0.7 


3.1xl0 5 y 


38 C1 


Spallation 


116.6 


4.9 


37.24 m 


35 Ar 


(n,6n') 


7.5 


6.0 


1.8 s 


37 Ar 


(n,4n') 


22.3 


0.8 


35 d 


39 Ar 


(n,2n') 


96.5 


0.57 


269 y 


41 Ar 


capture 


47.8 


2.5 


109 m 


38 K 


Spallation 


0.69 


5.9 


7.6 m 


40 R 


(P> n ) 


6.9 


1.3 


1.28 xlO 9 y 


Total 




1741 







a complicated simulation package. Analytic models are 
developed based on the physics processes that are then 
compared to the experimental results. Therefore, ana- 
lytic models can be used to compare and evaluate the 
results from a full Monte Carlo simulation. We develop 
analytic models below. 



1. Solid Angle 

When considering the interaction of particles in a de- 
tector, it is necessary to calculate the solid angle sub- 
tended a certain distance from the target. Using the 
work of H. Gotoh and H. Yagi |22j, the solid angle of 
a point particle was calculated at an arbitrary distance 
from the detector as well as the subtended solid angle 
from a fixed point throughout the detector material. For 
an assumed LBNE detector of dimensions (2w x 21 x h) 
= (45.6 x 22.4 x 14) m 3 , the solid angle from a point 

( x p: Yp: x p) IS 



fi = arctan 



arctan 



arctan 



- arctan 



(x p + w){y p + I) 



Zp((x p 4- w) 2 - 






- z 2 


1/2 


(Xp + 


w){y P 


-I) 






Zp((x p + w) 2 - 


+" (y P - 




- z 2 


1/2 


(Xp — 


w){y P 


+ 1) 






Zp{{x p - w) 2 - 


+ 




- z 2 


1/2 


(Xp — 


w)(y p 


-0 







B. Development of Analytic Models 



Zp({x p - w) 2 + (y p - I) 2 + z 2)Va 



(2) 



It is of general interest to have analytic models that 
can be used to estimate the cosmogenic production rates 
without doing a complicated simulation procedure. Even 
after a campaign of Monte Carlo simulation, it is neces- 
sary to assess whether a correct result was delivered by 



With this equation the average solid angle for a sheet of 
particles generated immediately above the detector was 
calculated to be 1.74 sr, and the differential muon flux 
from [23| could be used in the calculation of the total 
muon flux for a given level at Homestake mine. 
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2. Capture Rate for Negative Muons 

The stopping muons can be captured by an argon nu- 
cleus to generate radioactive isotopes such as 40 C1, etc. 
There are two sources of stopping muons for a given 
depth in an underground laboratory: 1) the through- 
going muons come to the end of their energy range and 
2) the secondary muons generated locally by the primary 
muons and their daughter pions. The capture rate for 
muons as a function of the depth of the detector can be 
calculated using the following equation: 



R capt 



ftp, ' fcapt ' fc ' fch ' fgi 

where the contributing terms are as follows: 



(3) 



1. is the stopping muon rate (derived below). 
2- fcapt is the fraction of muons that are captured, 

^lifetime 



with T capt 



fcapt 

Tcapt 

muon capture time in argon, 
11 1 

TO 



(4) 



(5) 



Tcapt ^lifetime 

Tufetime = 537±32 ns in argon [24[ and To = life- 
time of muons in a vacuum. Here, f cap t = 0.76. 

3. f c is the elemental fraction of the target in the com- 
pound 



fc 



Zi 



a i ■ Zi ' 

For the LBNE detector L = 1. 



(6) 



4. i c h is the charge ratio of negative muons to total 
muons, 



ch 



ii 



(7) 



It is 0.44 on the surface, but it is assumed that 
underground the fraction will be similar. 

5. f g is the fraction of reactions that occur with 40 CI 
in the ground state. Two fractions are considered: 
'gb = 0.2 (from the Geant4 



f ga = 0.0712 [23 and f, 



simulation) . 

The stopping muon rate can be calculated with the 
following formula: 



R-f 



scale 



(8) 



where R^ is the through-going muon rate, R is the ratio 
of stopping to through-going muons and, f scale, the scale 



factor 



' area lQQg-cm 



scales from a 1 m.w.e detector to a 



larger size and was calculated to be 19.5 for the simulated 
LBNE detector. 

The through-going muon rate is defined as 



Ri 



n-s a 



(9) 



A„ is the differential through- 



where Q is the solid angle, 

going muon flux, and S ar ea is the area through which the 
through-going muons transvers the detector. Using the 
measured differentiated through-going muon flux from 
F.E. Gray et al. [23| and the average solid angle, the 
calculated through-going muon rate is 4.2 xlO 6 per day 
at 800 ft and 5910 per day at 4850 ft. 

Two equations were analyzed for the ratio of stopping 
to through-going muons. The first, 



Ri = 



0.3 

<E fl > 



■5.7 -10" 



(10) 



was proposed by Chudakov et al. in Ref. [26| using ex- 
perimental data, and considers both the through-going 
portion as well as the muons produced by cascades with 



< E^> defined as 



<E p > -- 



t (l -exp(-bh)) 
7m "2 



(11) 



where no varies from 0.4-0.75 depending on depth (here, 
no = 0.4 was used to fit the equation best with simulated 
data), 7 M = 3.77 [H, e M = 693 GeV, b = 0.4 km.w.e, 
and h is the depth in km.w.e [30l |. The second equation 
considered, 



i? 2 = 



7/^ 



AE ■ exp(|) 
[esp(4) - % 



(12) 



is a parametrization from (2lJ, where 7 P = 3.77 [3(j, £ = 
2.5 km.w.e., AE w a Ax, a = 0.268 GeV/km.w.e. (for 
E M > 1000 GeV) and Ax = lOOg cm" 2 , h is the depth 
in km.w.e, and e M = 618 GeV [UH. 

At each depth, two rates were calculated for stopping 
muons and compared to the simulated values. These can 
be seen in Table [V] 



TABLE V: Comparison of calculated stopping muon rates to 
simulated results. 





Level 


Process 


800 


4850 


RjLi (p er da y) 

R^,2 (per day) 
Sim. (per day) 


3.06 xl0 b 
5.45xl0 5 
3.69 xlO 5 


251 
233 
173 



The stopping muon rates are within 30% of the sim- 
ulated values. The total capture rate of negative muons 
was calculated using two equations for the ratio of stop- 
ping to through-going muons (R x and R2) and the two 
different ground state fractions (f 90 and f g f,) resulting in 
four values (Table IVT]). 
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TABLE VI: Comparison of muon capture rates. 



exp(- 



<L> 



)], is the probability that the remaining neu- 





Level 


Process 


800 


4850 


r mX (p er da y) 
KTa (p er da y) 


7286 
1.30 xlO 4 


5.97 
5.55 


R ^i6 (p er da y) 
(p er da y) 


2.05 xlO 4 
3.64xl0 4 


16.8 
15.6 


Sim. (per day) 


2.73 xlO 4 


17.5 



The capture rates and simulated values were plotted 
as a function of depth (km.w.e) in Fig. [5] using the in- 
tegrated neutron flux (Eq.(13) Ref. [13]) instead of the 
differential flux at the specific levels. As it is apparent 
in Fig. El using R™fj to calculate the muon capture rate 
has the best agreement to Monte Carlo. 




5 6 7 

Depth (km.w.e.) 



FIG. 8: Capture rates as a function of depth. 



3. Production Rate via (n,p) Reaction 



trons will undergo (n,p) reaction and produce the 40 C1 
background. The remaining terms, p .™ L> , are character- 
istics of the detector. 

In calculating the average path length, the angular 
dependence of the neutrons (and their parent particle 
muons) was considered using the methodology of Jostlein 
and McDonald 31]. A -^-g dependence was used in the 
calculation. 

The mean free path was calculated by 



(14) 



P ■ N a ■ a 



where A is the mass number of the target nucleus, p 
is the density, N a is Avogadro's number, and a is the 
cross section. For the calculation, the flux weighted cross 
section, a, was used 



(15) 



where flux and cross section of the same energy, i, are 
summed together. The cross sections from Geant4 were 
used for this calculation. It is important to use the flux 
weighted cross sections to calculate the mean free path 
in order to account for the entire spectrum of possible 
neutron energies. Unfortunately the Geant4 cross section 
data does not exceed 100 MeV, however, the calculated 
values should be accurate to within a factor of two. 

Using the above formulas, we have calculated the av- 
erage path length, the flux weighted cross sections, and 
the mean free path for neutrons at the levels of 800 ft 
and 4850 ft. Table IVTA shows the results. As can be seen 
in Table I VIII the values are similar as anticipated. 



Another important reaction that contributes to the 
production of 40 CI in the detector is the (n,p) reaction, 
40 Ar(n,p) 40 Cl. To calculate the production rate the fol- 
lowing formula can be used: 



$„(£„) • exp{- 



<L> 

^total 



)[1 - exp(-f±>-)] ■ m 



(",p> 



p- < L> 



(13) 

where m is the mass of the detector, p is the density of 
the detector medium, < L > is the average path length, 
&n(E n ) is the integrated neutron flux, Xtotai is the mean 
free path considering all neutron disappearance reaction 
channels except (n,p), and X( n ,p) is the mean free path 
of the (n,p) reaction channel. 

The first component of the (n,p) production rate, 
$ n (E n ) • exp(— ^ L> ), is the probability that neutrons 



will survive all reactions except the (n,p) reaction, re- 
ducing the flux and availability of neutrons for the 
production of 40 C1. Similarly, the second term, [1 — 



TABLE VII: The calculated parameters for Eq. [T3J 





Path length 


Cross section 


Mean free path 


Level 


< L > (cm) 


ff(n,p) (b) 


Ototal (b) 


X(n,p) ( cm ) 


Xtotai (cm) 


800 


572.2 


0.0209 


0.181 


2280 


262.9 


4850 


570.8 


0.0209 


0.181 


2280 


262.9 



Two main sources of neutrons are considered for both 
800 ft and 4850 ft: 1) the muon-induced neutrons en- 
tering the detector from the experimental hall (Source I 
neutrons) and 2) the muon-induced neutrons in the de- 
tector when muons pass through the target (Source II 
neutrons) . The final calculated (n,p) production rate for 
both levels are listed in Table IVTlTl with the relevant neu- 
tron flux in the detector above the (n,p) reaction thresh- 
old. 

As shown in Table [VlTTl the predicted production rates 
using the analytic models agree well with the Geant4 
simulation. 
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TABLE VIII: The calculated (n,p) production rates. 







^(-CJn) (cm S ) 


4U C1 rate per day 


Source I 
Neutrons 


Level 
800 
4850 


Flux 
8.8xl0" 7 
1.66xl0~ 10 


Analytic models 
4784 
9.07 


Geant4 
3667 
9.3 


Source II 
Neutrons 


Level 
800 
4850 


Flux 
1.4xl0" 6 
6.0xl0~ 10 


Analytic models 
66978 
54 


Geant4 
40587 
54 



C. Scaling Function 

Muon-induced processes and the cosmogenic radioac- 
tivity production depends strongly on the target and it 
must be evaluated individually for the experiment. How- 
ever, the production rate is proportional to muon flux, 
or neutron flux, and their interaction cross-section. The 
energy dependence of the total cross-section for all muon- 
induced radio-isotopes in the scintillator was evaluated 
assuming the power law [32j 

otrf^) a ££, (16) 

where a varies from 0.50 to 0.93 with a weight mean value 
< a > = 0.73±0.10 [33j. For a given target atoms N 
and the cross-section ao at the Earth surface where the 
average muon energy is about 4 GeV, the muon-induced 
cosmogenic radioactivity (Rj SO ) depends on the differen- 
tial muon energy spectrum dN^/dE^ at the experimental 
site at a depth ho, 

^-^f'lf^i^' < 17 > 

As a simplification, the production rate is written as a 
function of the average muon energy < > at a depth 
ho M- 

Rise = a N<T 4 Oev C^^y )"^, (18) 

where </> M is the total muon flux at the experimental site 
and 0o.73 = 0.87 ± 0.03 is the correction factor for the 
averaging of E M [33| ■ For a given detector target and a 
depth, the cosmogenic production rate as a function of 
depth is thus obtained 

R iso (unknown) _ E^ unkn own , a j ^{unknown) 
R iso (known) < E^ known > ^(known) 

(19) 

III. DEPTH-SENSITIVITY RELATION 

LBNE is an extremely rich physics program that will 
measure neutrino properties using a neutrino beam. In 
addition, LBNE will also measure supernova neutrinos 
and proton decays. Each of these physics channels has 
its unique signal region in terms of energy distribution. 



For example, measuring the parameters of neutrino oscil- 
lation with a neutrino beam has an energy region from 1 
to 8 GeV while the energy region of proton decay ranges 
from 100 to 938 MeV. The signal of supernova neutrinos 
resides between 5 to 50 MeV. It is difficult to establish a 
depth-sensitivity relation for all physics channels using a 
single plot. We elaborate the muon-induced backgrounds 
for each of physics channels below. 

A. Muon and the Muon-induced Neutron Rates 

Utilizing the above formulas, the muon and the muon- 
induced neutron rates as a function of depth can be 
shown in Fig. [9l It is clear that the event rates de- 
crease rapidly when the depth increases. The fluctuation 
of muon and the muon-induced neutron rates from sea- 
sonal variation can result in backgrounds for all physics 
channels depending on the depth. 




Depth (km.w.e.) 



FIG. 9: Muon and muon-induced neutron rates as a function 
of depth. Shown is for the energy deposition greater than 5 
MeV in the detector. 



B. Cosmogenic Production Rates 

The cosmogenic production rate also decreases with 
increasing of depth. As an example, Fig. [TU] shows the 
40 C1 production rate as a function of depth. The cosmo- 
genic produced radioactive isotopes can be backgrounds 
for the detection of relic supernova neutrinos. 

C. 7To Production Rates 

The neutrino beam-induced natural current and charge 
current 7To productions are important backgrounds to 
the v e appearance. However, ttq can also be produced 
by muons and muon-induced neutrons in the detector. 
Fig. [TT] shows the energy deposition from 7To events cre- 
ated by fast neutrons in the detector. It is worth mention- 
ing that there are also multiple ttq events along the neu- 
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Induced by neutrons 
Induced by survivial muons 




6 7 
Depth (km.w.e.) 
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CO 
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Muon-induced n Q 
Fast neutron induced it. 




6 7 
Depth (km.w.e.) 



FIG. 10: CI production rates produced by (n,p) reaction as 
a function of depth. 



FIG. 12: 7to production rates as a function of depth. 



tron track. The production of ttq as a function of depth is 




10' 10" 
Energy Deposition (MeV) 



FIG. 11: 7ro events in the detector at the 4850-ft level. 

shown in Fig. [12J It is clear that the ttq produced by fast 
neutrons entering the detector can be a significant back- 
ground to the beam physics program. This is because 
the fluctuation of the n production due to the seasonal 
variation of muon flux and neutron flux results in a level 
of 22 tto events per year in the detector. In addition, the 
statistical fluctuation of the ttq has a similar level of 33 
7To events. Adding both in quadratic, the level of fluc- 
tuation in the production of ttq can be about 40 events 
per year. This is significant even at the 4850-ft level. 
Depending on the capability of discriminating single tt 
events between the beam neutrino-induced and the cos- 
mic neutron-induced, a greater depth (7400-ft level) can 
be an effective option to further reduce this background. 



D. Backgrounds for v e Appearance 

The v e appearance is essential to the measurements of 
neutrino properties with a neutrino beam. The antici- 
pated signal in a 20 kton LAr detector is about 75 events 
per year [2l| . The muon-induced backgrounds can be 



measured with the beam-off. However, the fluctuation 
of background events resulting from statistical and sea- 
sonal variation is a main source of background. This 
background as a function of depth is shown in Fig. 1131 
As can be seen in Fig. [TSJ the depth must be greater than 
4.0 km.w.e. in order to have a reasonable measurement 
of u e appearance at SURF. 



Fast neutron induced background 

v„ Appearance 

All Other Backgrounds 




6 7 
Depth (km.w.e.) 



FIG. 13: Muon-induced background as a function of depth. 
Shown is the average value of statistical fluctuation induced 
by fast neutrons. Note that signal and the other total back- 
grounds were estimated using Ref. [2lj . 



IV. CONCLUSION 

We evaluate the muon-induced background as a func- 
tion of depth for a long baseline neutrino experiment with 
liquid argon as the target at the SURF at Homestakc. 
Both Geant4 simulations and analytic methods are em- 
ployed in the evaluation of background event rate in the 
region of interest using the available muon and neutron 
energy spectra from Ref. [l7| and the measured muon 
flux from Ref. [23[ . The muon and muon-induced neutron 
rates are calculated for the energy greater than 5 MeV. 
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The cosmogenic processes are discussed in detail for nega- 
tive muon capture, (n,p), and (p,n) reactions. The dom- 
inant backgrounds are from both stopping muons and 
the high energy neutrons generated by muons in the sur- 
rounding materials and in the target. In summary, the 
background sources are: 1) the negative stopping muon 
capture; 2) the neutron (n,p), (n,d),(n,t), (n,a), etc; and 
3) the muon-induced showers. The cosmogenic produc- 
tion rate as a function of depth is evaluated for negative 
stopping muon capture and fast neutrons separately. As 
can be seen from the above discussion, the cosmogenic 
production rate induced by the muon-induced processes 
reduces significantly when the depth is larger than 4.0 
km.w.e. It is clear from Fig. [8] and Fig. QUI the 40 C1 pro- 
duction is less than 100 per day when the depth is greater 
than 4.0 km.w.e. We conclude the following: 

• The 800-ft level is too shallow for a beam-related 
neutrino physics because the muon rate (^88 Hz) 
is too high. Though LAr TPC is better in iden- 
tifying 7To events comparing to a water Cerenkov 
detector [34|], the ttq produced by muons seen in 
Fig.rj2]can be a significant background. 

• The 800-ft level cannot be used to detect supernova 
neutrinos for a galactic supernova neutrino burst 
with a time window of 30 seconds at 10 kpc. This 
is because the expected charge current events in a 
20 kton detector used in the simulation from such a 
burst are about 1300 events 20]. They corresponds 
to about 44 Hz. However, the total muon-induced 
rate is about 88 Hz at this level. The signal is im- 
mersed in background. In addition, the frequency 
of galactic supernova occurs at a level of once per 
50 years (35| . 

• The muon-induced processes can be significant 
background for an argon-based detector in the de- 



tection of relic supernova neutrinos with a depth 
less than 4.0 km.w.e. With a 20 kton detector, we 
expect less than 30 events per year from relic super- 
nova neutrinos 36] . The relic supernova neutrinos 
can be detected with the accumulation of the de- 
tector lifetime. 



• Fig. [T3] shows a depth requirement for the v e ap- 
pearance from a neutrino beam. It is clear that a 
meaningful measurement of CP violation can only 
be accomplished when the depth is greater than 4.0 
km.w.e. 



• Finally, positioning a LAr detector near surface, i.e. 
at NOVA depths, will increase the backgrounds by 
three orders of magnitude compromising the direct 
v program. 

Therefore, a depth larger than 4.0 km.w.e. is needed 
for an argon-based detector. The 4850-ft level at Home- 
stake mine can be an appropriate home for this detector. 
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